Scanning Tunneling Spectroscopy of the superconducting proximity effect in a diluted 

ferromagnetic alloy 
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We studied the proximity effect between a superconductor (Nb) and a diluted ferromagnetic alloy 
(CuNi) in a bilayer geometry. We measured the local density of states on top of the ferromagnetic 
layer, which thickness varies on each sample, with a very low temperature Scanning Tunneling Mi- 
croscope. The measured spectra display a very high homogeneity. The analysis of the experimental 
data shows the need to take into account an additional scattering mechanism. By including in 
the Usadel equations the effect of the spin relaxation in the ferromagnetic alloy, we obtain a good 
description of the experimental data. 



PACS numbers: 74.45.+C, 74.78.Fk 



I. INTRODUCTION 



Hybrid systems made of a ferromagnetic metal (F) in 
contact with a superconductor (S) exhibit a rich vari- 
ety of physical effects. At the interface, the Andreev 
reflection 0, Q is responsible for the appearance of the 
proximity superconductivity in the ferromagnetic metal. 
In the ferromagnetic metal, the reflection of a spin up 
(down) electron into a spin up (down) hole correlates 
the occupancy of the two electron states of opposite spin 
bands. This process is equivalent to the creation of an 
electron pair (an Andreev pair). In the framework of a 
Stoner model for the ferromagnetic metal, spin up and 
spin down energy bands are shifted by the exchange field 
E ex . This implies that the paired electrons states have a 
different wave- vector. At the Fermi level, the wave- vector 
mismatch is 8k — kpE ex /Ep where kp and Ep are the 
Fermi wave-vector and energy. The Andreev pairs have 
therefore a non-zero momentum, which induces an oscil- 
lation of the pair density as a function of the distance 
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to the interface 0. [J Q In a diffusive ferromagnetic 
metal, the oscillation goes together with a decay on the 
same length scale, i.e. the ferromagnetic coherence length 
£p = \JhDpjE ex where Dp is the diffusion coefficient in 
F. The electron energy spreading contributes also to the 
wave-vector mismatch, but usually in a negligible way, 
since the exchange field E ex is usually much larger than 
the thermal energy ksT . 

Depending on the distance to the interface, the super- 
conducting wave-function phase changes from to n. In 
a S/F/S junction, a 0-phase or a 7r-phase regime is ex- 
pected, depending on the ferromagnetic metal length as 
compared to the coherence length £p . The crossover be- 
tween a 0-phase and a 7r-phase regime was observed in 
Nb/CuNi/Nb junctions laQ- There, the thermal energy 
is not negligible as compared to the exchange field and a 
sign change of the supercurrent was observed when the 
temperature was varied. In Nb/PdNi/Nb junctions, the 
critical current does not change sign with temperature 
but still exhibits an oscillating behavior as a function of 
the ferromagnet thickness A 7r-regime was also de- 
tected in an interference device made of five 7r-j unctions 
and in a SQUID made of one 0-j unction and one ir- 
junction |T(ij]. 

The 7r-regime can also be observed through density 
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FIG. 1: Spectrum measured on a plain Nb layer (full line) and 
a fit with a BCS calculation (dashed line). The fit parameters 
are the energy gap A = 1.5 meV and an effective electronic 
temperature of f 70 mK. 

of states measurements in a S/F junction. The first 
observation was performed with solid tunnel junctions 
deposited on top of Nb/PdNi bilayers ^lj. For a cer- 
tain range of PdNi thickness, the density of states ap- 
peared as reversed, with a maximum at the Fermi level 
and a minimum at the gap energy. The data were fit- 
ted with the quasi-classical theory ^^|- The density of 
states modulation appeared to be surprisingly small, and 
shifted to large ferromagnetic alloy thickness. The lat- 
ter fact was accounted for by the possible existence of 
a magnetically-dead layer at the interface. The spatial 
homogeneity of the observed density of states was not 
investigated because of the intrinsic spatial averaging of 
the technique. This leaves open questions that motivated 
us to study similar samples with a very low tempera- 
ture scanning tunneling microscope (STM). The spatial 
resolution should constitute a clear benefit of this new 
technique, while the energy resolution is comparable. 



II. EXPERIMENTAL STUDY 

Our study relies on six different wafers made with the 
same schematics. On every silicon wafer, we deposited 
first a thick layer of Nb and then a thin layer of a CuNi 
alloy. The CuNi target composition was CU50M50. The 
sputtering conditions were identical to that of Ref. 0, 
where Nb/CuNi/Nb junctions critical current was stud- 
ied. The Nb layer thickness of 200 nm (at the wafer 
center) was chosen to be much larger than the Ginzburg- 
Landau coherence length which is estimated to be 12 nm. 
The CuNi film thickness varies on each sample (from 2.5 
to 20 nm at the wafer center). Because of the charac- 
teristics of the chamber, there is a CuNi film thickness 
gradient from the center to the sides of the wafer. This 
was taken into account in every indicated thickness in 
the following. 



FIG. 2: Superposition of 20 tunneling spectra measured on a 
sample with a 3.3 nm-thick CuNi film. Spectra were acquired 
every 11 nm along a single line. The tunnel conductance is 
of the order of 100 nS. Left insert : 400 x 400 nm 2 image of 
a sample with a CuNi film thickness of 3.3 nm. The black to 
white scale is lnm. Right insert : Schematics of the measure- 
ment geometry. 



CuNi (like PdNi) is a diluted ferromagnetic alloy with a 
relatively small exchange field. This brings the advantage 
of a sizable coherence length £p of the order of several nm. 
On the other hand, the ferromagnetic properties depend 
strongly on the exact composition of the alloy. The Cu 
and Ni concentration was measured on one sample by 
Rutherford BackScattering to be CU52M48. The Curie 
temperature of similar films was measured to be about 20 
K and the saturation magnetization was extrapolated to 
5.10 A.m~ , which corresponds to a magnetic moment 
of 0.06 \XB-atom~ 1 (see Ref. for more details). As 
both the Curie temperature and the magnetic moment 
depend linearly on the concentration, an extrapolation 
[l3| for 48% indicates a value E ex = 12 meV. Taking 
into account a mean free path in CuNi of about 1 nm, 
which corresponds to a electron diffusion coefficient Dp 
of 5 cm 2 /s, this gives an estimated coherence length £p 
= 5.0 nm. 

An important property of our CuNi films is the ab- 
sence of hysteresis or remanent magnetization f?). This 
property indicates a very weak ferromagnctism close to 
the superparamagnetic behavior that indeed occurs at 
45% of Ni in bulk alloys ^4[. In this regime, the mag- 
netic correlations are short range and the direction of the 
magnetization varies spatially without domain wall. The 
characteristic length of these magnetic fluctuations is a 
few hundred atoms, close to the superconducting coher- 
ence length, so that the exchange energy should appear 
as uniform in terms of induced superconductivity. 

We used a very low temperature (60 mK) Scanning 
Tunneling Microscope (STM) developed in our labora- 
tory |15|. Local spectroscopy measurements were per- 
formed at a temperature below 100 mK and with a typi- 
cal tunnel resistance value of 10 to 25 Mfl. The measured 
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FIG. 3: Tunneling spectra (full lines) measured at the surface 
of three Nb/CuNi bilayers with different CuNi thicknesses 
2.5 nm (circles), 3.3 nm (squares) and 5.3 nm (diamonds), 
compared with calculated tunneling spectra (dotted lines). 
The calculation parameters are £f = 15 nm, T = 170 mK 
and respectively L s = 3.4 ; 2.6 and 3.0 nm, A = 1.5 ; 1.45 
and 1.4 meV. 



I(V) characteristics (I is the tunnel current, V the volt- 
age applied on the tip) were numerically differentiated 
so that we obtain the differential conductance dl/dV. In 
the limit of a zero temperature, this quantity is directly 
proportional to the local density of states (LDOS) of the 
sample. In order to test the energy resolution of the 
STM, we measured tunneling spectra on a plain Nb film, 
see Fig. Q] The spectra is fitted with a thermally smeared 
BCS density of states, without using any inelastic scat- 
tering parameter. The fit parameters are the energy gap 
A = 1.50 meV and an (effective) electron temperature T 
= 170 mK. This value represents an improvement as com- 
pared to our earlier experiments |lj|. Nevertheless, the 
difference with the measured sample temperature (below 
100 mK) means that some non-thermal noise is present 
in the STM. 

Fig. |21 shows the superposition of 20 spectra all taken 
on the 3.3 nm of CuNi sample while scanning along a 
line with a step size of 11 nm. All spectra are almost 
perfectly superimposed. A similar behavior has been ob- 
served on every sample. This means that there are no 
inhomogeneities in the LDOS up to a scale much larger 
than £f. This is possible because of the excellent rough- 
ness of the samples, of the order 0.5 nm, and the good 
stability of the CuNi alloy surface. Fig. [2]shows a typical 
image, of size 400 x 400 nm 2 . No clear feature can be 
observed. 

We measured about ten different samples, some of 
them were taken from the same wafer but at a differ- 
ent place. A selection of data is shown in Fig. [3]for three 
different thicknesses x of CuNi. The comparison of these 
spectra with the Fig. ^Nb data demonstrates that the 




FIG. 4: Normalized differential conductance vs. bias voltage 
on top of a Nb (100 nm)/ CuNi (5.1 nm) bilayer. Note the 
scale in dl/dV. Experimental noise can be observed. Two 
dashed lines are drawn at V = ±A/e for clarity. 



smearing of the spectra is not due to a lack of experi- 
mental resolution but is a physical effect. For a small 
thickness of 2.5 nm, there is a sharp rise to a peak in 
the measured LDOS at V — ±A/e. One can note that 
A ~ 1.5 meV, which is very close to the value for bulk 
Nb (1.6 meV). With a CuNi film thickness of 3.3 nm, 
the density of states minimum at the Fermi level is less 
pronounced while the maximum position does not change 
significantly. This trend appears more strongly in the 5.3 
nm sample. As will be discussed below, the behavior of 
the measured tunneling spectra is not fully monotonous 
as a function of the CuNi film thickness. For a CuNi 
thickness larger than 6 nm, where the 0-tt crossover is 
expected, the density of states is flat within our exper- 
imental accuracy of about 1 %. This limitation comes 
from the imperfect junction stability that restricts the 
time averaging possible in a STM experiment. 

In one experiment with a CuNi thickness of 5.1 nm, we 
observed a reversed density of states (see Fig. 0J. This 
kind of spectra is reminiscent of the 7r regime, which is 
predicted to appear in a S-F junction at a F metal thick- 
ness x > 7t£f/2. The background is not as flat as usual. 
A clear structure is observed at voltages close to ±A = 
1.5 meV (dashed lines), which confirms that this is actu- 
ally an effect of the proximity with the superconductor. 
This behavior was reproduced in more than ten succes- 
sive spectra acquired at different positions, but eventu- 
ally disappeared as the tip was scanned further. We be- 
lieve that these spectra did rely on an uncontrolled sam- 
ple configuration below the tip. For instance the tip may 
have picked up some CuNi cluster at the top of the sam- 
ple, so that we performed spectroscopy through a S-F- 
I-F junction (I for "insulator"). Additional experiments 
in other samples within the same thickness range did not 
reproduce the same behavior. 



III. THEORY 

Let us now discuss the theoretical description of the su- 
perconducting proximity effect in a ferromagnetic metal. 
We will consider the diffusive regime where the elastic 
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mean free path is supposed to be smaller than every 
characteristic length scale, including the ferromagnetic 
metal length. In the diffusive regime, the quasi-classical 
theory for the inhomogenous superconductors reduces to 
the Usadel equation. Using the usual parametrization of 
the normal and anomalous Green functions G — cos 
and F = sin 8, the Usadel equation is written for u> > 
as: 



Dp d 2 e 



+ (u> + iE ex ) sinG + — sin 9 cos 8 = 0. (1) 

OX 1 T s 



eters us = u>/E ex , t s = T s E ex and x = x/£,f, Eq. ^ rea( is: 



1 <9 2 8 „ 1 

+ (u + i)sm& + — sin8cos8 = 0. (2) 

2 ox z T s 



For an infinite ferromagnet length, the integration of 
(J2J using the boundary conditions: (x — > oo) = and 

-^-zr ) = allows us to find the exact expression : 
ox / - _ 



Here we considered a case where the physical quantities 
depend only on the coordinate perpendicular to the lay- 
ers. Compared to the usual expression of the Usadel 
equation, the magnetic scattering in this case is incorpo- 
rated by replacing u by ui + G/t s |l6j . where t s is the 
magnetic scattering time. Note that we assume the pres- 
ence of a relatively strong magnetic uniaxial anisotropy 
in the ferromagnetic layer. This permits us to neglect 
the magnetic scattering in the plane perpendicular to 
the easy axis (which mixes up and down spins in the 
Green functions). Introducing the dimensionless param- 



G = cosG = 1 



where k 2 = 
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and g is defined as: 
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Go = G (x — 0) is the normal Green function in the S 
layer, that is 



Gq — 



VA2 



- UJ* 



(6) 



This means that we neglected the inverse proximity effect 
and assumed that the pair amplitude at the interface is 
equal to the one in the bulk. The density of state is 
deduced from G performing the analytical continuation: 



and the decay length are given by the coherence length 
7r£p. As can be seen from Eq. the spin relaxation 
influences the oscillations of the LDOS in two ways: it 
diminishes the damping characteristic length of the oscil- 
lations and it also increases the period of the oscillations. 
In the regime f s C 1, the period is multiplied by the fac- 
tor (1+2^) which is larger than one. The decay length is 
divided by the same factor. In the opposite regime, the 
period is ■k\[2^ 1 f j L s while the decay length is the spin 
relaxation length L s = yJDsr s . 



N (e) = N(0)ReG(x,u -> ie) 



(7) 



where e = E — Ep is the energy referenced to the Fermi 
level. With no spin relaxation, both the oscillation period 



For a finite ferromagnet length dp, the boundary con- 
ditions become Q (x — dp ) = Qd F and ( -— I = 0. 

^ ' \ox J 

A first integral of (J2J gives the following equation: 



J 



dx = 



de 



2 (cos Q dF - cos 8) (2 (u + i) + i (cos Q dp + cos 6) 



(8) 



where Qd F is the parametrization function at dp. In the case of Qd F 1, the integration of © can be analytically 
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performed and we obtain: 



: exp 



i + — + ujcIf 



l-p 2 sin 2 (^) 



cos(^) 



\ Jl-p 2 sin 2 (^)+cos(%) 



(9) 



where p 2 = — . The normal Green function at the 

1 + IT S 

end of the F layer is deduced from and then, the 
LDOS can be determined with the analytical continua- 
tion. 



IV. DISCUSSION 

Fig. El shows the LDOS change at the Fermi level 77 as 
a function of the samples CuNi film thickness. A value 
of 100 % for 77 would mean a zero density of states at the 
Fermi level, while the normal state gives 77 = 0. There 
is a significant scatter in the data, which cannot be ex- 
plained by the error bars on the CuNi thickness. The 
overall behavior is a sharp decrease on a length scale be- 
low 3 nm, with no visible 7r-regime. This decay length is 
significantly smaller than the expected coherence length 
This shows that our experimental results do not fol- 
low the predictions of the usual theory for F-S junctions. 

We added with a different symbol the data point of the 
7r-regime-like spectra observation. The reversed density 
of states observed in one experiment may be understood 
as a special realization of a 7r-regime due to the tun- 
neling through a nanometric ferromagnetic grain weakly 
coupled to the sample surface. It was actually predicted 
that a weakly transparent interface in a S-I-F junction 
can induce a 7r-regime in a very thin ferromagnetic layer 
[r?| . The physical explanation is that the pair amplitude 
drop at the interface mimics the effect of a ferromagnetic 
metal of some thickness. No spin relaxation is expected 
in the barrier, so that the 7r-regime could be of stronger 
amplitude than in the transparent interface case. 

When the F layer is rather small (2.5 nm), 77 is higher 
than 70%. If one extrapolates 77 at the interface (Onm 
of CuNi), one can note that it would reach a value close 
to 100%. This implies that the interface transparency is 
very good, since a moderate interface transparency would 
result in a pair amplitude mismatch, and consequently a 
density of states aperture mismatch. 

We interpret our data by taking into account a large 
spin relaxation in the CuNi film. This effect was invoked 
for fitting the critical current measurements results in 
S/F/S junctions based on the same CuNi layers 0. In 
qualitative terms, a spin relaxation length L s shorter 
than the coherence length £p will be the leading term 
in the decay of the pair amplitude with the distance to 
the interface. The spin relaxation length estimated in 
L s = 2.7 nm is at least in qualitative agreement with our 
experimental findings. 
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FIG. 5: Evolution of the change of Fermi level density of 
states as a function of the samples CuNi film thickness. A 
value of 100% for the parameter 77 would correspond to a gap 
in the density of states. Data are shown with a filled symbol, 
except for the reversed density of states data (see text) shown 
with an open symbol. 

Figure shows an experimental spectra of the 4.1 nm 
sample together with several calculated spectra. The 
density of states was calculated using Formula|§land con- 
voluted with a thermal window in order to obtain the 
predicted tunneling spectra. The energy gap was taken 
as A = 1.37 meV and the temperature taken equal to 
170 mK. The two free parameters were the ferromagnetic 
metal coherence length £p and the spin relaxation length 
L s . If we take the parameters values that were derived 
from Ref. critical current measurements (L s =2.5 nm 
and £p=4.1 nm), one obtains a calculated spectra (trian- 
gle symbols) with an modulation amplitude that is much 
smaller than in the experimental spectra. Increasing the 
spin relaxation length (L s = 4.5 nm, cross symbols) en- 
hances the calculated amplitude but the agreement with 
the experimental data remains very poor. 

In contrast, a smaller exchange field gives a good fit 
to the experimental data. We obtain the best agreement 
with a spin relaxation length L s = 3.6 nm and a co- 
herence length £p = 15 nm (circle symbols). The latter 
value corresponds to an exchange field E ex of 1.5 meV, 
well below the expected value. We interpret this as a 
signature of the weakened magnetism of the CuNi mate- 
rial in the vicinity of the Nb interface. Since we inves- 
tigated mostly very thin CuNi films (below than 6 nm), 
this statement is compatible with the observation of a 7r 
regime crossover at 17 nm 0, driven by the larger ex- 
change field of thicker films. A magnetically-dead layer 
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FIG. 6: Comparison of calculated density of states spectra 
(thin lines) with the experimental for the 4.1 nm sample (thick 
line). The calculation parameters are : A = 1.37 meV, T = 
170 mK, £f ; L a ) = (15 ; 5) in nm with square symbols ; (15 
; 3.5) with circles ; (4.1 ; 2.5) with triangles and (4.1 ; 4.5) 
with diamonds. 

of 1.5 nm at the interface was invoked in PdNi in or- 
der to interpret quantitatively density of states measure- 
ments 11]. In the present case of CuNi, we claim that 
the magnetism of the diluted ferromagnetic alloy is not 
destroyed but only reduced in the vicinity of the Nb inter- 
face. The spin relaxation length L s = 3.6 nm extracted 
from this fitting procedure is in fair agreement with Rcf. 
and also 0. Let us note that the calculated spec- 
tra depend very strongly on the value of spin relaxation 
length. For comparison, Fig. Elalso shows the calculated 
spectra with a slightly larger spin relaxation length L s 
— 5 nm (square symbols). The modulation amplitude 
is strongly enhanced and the agreement with the data 
becomes much poorer. 

Fig. |3 shows the comparison between experimental 
data and the related best fit calculated curves. The co- 
herence length parameter was kept constant with £p = 15 
nm, while the spin relaxation length and the energy gap 
were adjusted in every case. The agreement is very good. 
The 2.5 nm data is less accurately described, maybe be- 
cause a lack of validity of the diffusive regime approxi- 
mation. In this sample, the CuNi thickness is actually 
of the order of the elastic mean free path. A small vari- 
ation of about ± 12 % in the spin relaxation length is 
sufficient to account for the scattering of the data. This 
variation is present even between samples from the same 
wafer (for instance the 3.3 and 4.1 nm sample). It may be 



related to small differences in the chemical composition 
of the CuNi alloy (below the RBS sensitivity of 1%) on 
a scale larger than the scanning area of the STM. Such 
small differences can have a dramatic effect on the spin 
relaxation, especially in the concentration region close to 
the super-paramagnetic regime. 

The energy gap value extracted from the fitting proce- 
dure shows a tendency to decrease as the ferromagnetic 
metal length increases. This small effect may be under- 
stood as a small reduction of the energy gap at the F-S 
interface due to the inverse proximity effect. 



V. CONCLUSION 

In summary, we investigated experimentally the lo- 
cal density of states in a thin CuNi film in proximity 
with a Nb film. Our results are consistent with previ- 
ous studies and bring additional insight into this system. 
The measured tunneling spectra were very homogenous 
over the accessible sample area. This may be related to 
the fact that no magnetic domain structure is expected. 
The steady evolution of the measured spectra with the 
CuNi film thickness demonstrates that the interface is 
very transparent. 

The inclusion of the spin relaxation into the Usadel 
equation modifies the predicted LDOS in a significant 
way. Comparing our experimental data with theoretical 
prediction, we find that the spin relaxation has an essen- 
tial effect on the proximity superconductivity induced in 
a diluted ferromagnetic alloy like CuNi. This mechanism 
is responsible for the very small amplitude of the den- 
sity of states modulation in the 7r regime. The observed 
data scattering is accounted for by small sample to sam- 
ple variation of the spin relaxation strength. Moreover, 
we confirm that, in the vicinity of the Nb interface, the 
CuNi ferromagnetism is significantically reduced. 
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